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DE LA (N&FEGTUBE 




PROCEDE ET DISPOSITIF DE DETERMINATION DE LA POSITION 
RELATIVE DE DEUX POINTS, A BASE DE SIGNAUX DE 
POSITIONNEMENT PAR SATELLITES 



L'invention concerne la determination precise de la position 
relative de deux points qui peuvent etre eloignes de plusieurs dizaines de 
5 kilometres, a partir de signaux de positionnement par satellites. 

Le domaine (^application est celui des techniques necessitant de 
connaltre en trois dimensions et avec une precision centimetrique, la position 
d'un mobile par rapport a la position connue d'une station de reference 
eloignee de plusieurs dizaines de kilometres. On peut citer la geodesie la 
0 topographie, I'hydrographie, .... 

Pour la determination de la position relative d'un mobile par 
rapport a une station de reference, on utilise couramment des moyens de 
mesure de position par satellites, utilisant par exemple ies signaux radio emis 
par Ies satellites du systeme GPS (Global Positioning System) ou d'autres 
5 systemes analogues (systeme GLONASS, futur systeme GALILEO). 

Dans le systeme GPS, le signal emis par un satellite est code et 
on utilise le temps mis par ce signal pour atteindre le point a localiser pour 
determiner la distance entre ce satellite et ce point, de preference appelee 
pseudo-distance pour tenir compte d'erreurs de synchronisation entre 
I'horloge du satellite et celle de la station. Ces erreurs de synchronisation 
sont classiquement eliminees par calcul des lors qu'on recoit Ies signaux 
d'au moins quatre satellites differents. La determination de la distance entre 
le point a localiser et plusieurs satellites permet, connaissant Ies 
coordonnees geographiques des satellites, de calculer Ies coordonnees du 
point a localiser, le plus souvent des coordonnees exprimees en latitude, 
longitude et altitude dans un repere terrestre fixe. 

Pour determiner la position relative d'un mobile par rapport a une 
station de reference, on utilise une methode dite « GPS differentiel » qui 
consiste a localiser un point par rapport a une station de reference et non par 
rapport a un repere terrestre independant : en mettant en place un recepteur 
a la station de reference, on peut determiner, en utilisant Ies mesures faites a 
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la station et au mobile, la position relative du mobile par rapport a la station 
de reference. 

L'avantage de cette methode est qu'elle permet daccroitre la 
precision du positionnement. En effet les distorsions des mesures liees aux 
aleas de la propagation des signaux radiosatellitaires sont le plus souvent 
fortement correlees dans I'espace, et la connaissance precise de la position 
de la station de reference permet, en comparant les mesures faites a la 
station aux distances theoriques, de les compenser en grande partie. 

Le temps de propagation est determine d'une part en reference a 
un instant de repere du code pseudo-aleatoire qui module une frequence 
porteuse emise par le satellite, cet instant de repere du code permettant 
notamment de determiner la position approximative du mobile c'est-a-dire 
avec une precision de quelques metres a quelques dizaines de metres ; le 
5 temps de propagation est determine d'autre part en reference a la phase de 
la porteuse recue, la mesure sur la phase, moins bruitee que la mesure sur 
le code, permettant de determiner plus precisement la position du mobile 
c'est-a-dire avec une precision centimetrique a la condition toutefois qu'on 
leve Tambiguite sur le nombre de tours de phase puisque la phase ne peut 
>o etre a priori connue qu'a 2% pres, 2* correspondant a une distance egale a la 
longueur d'onde du signal radiofrequence emis par les satellites. 

Dans ce qui suit on ne s'interessera qu'aux mesures de phase, les 
mesures de position de code pouvant etre faites d'une maniere classique. On 
considerera done que les pseudo-distances fournies par le recepteur GPS du 
25 mobile ou de la station de reference sont essentiellement des valeurs 
numeriques de phase, une valeur de phase etant directement convertie en 
valeur de distance, connaissant la longueur d'onde du signal radio emis par 
les satellites. 

Le point central des techniques de positionnement centimetrique 
30 utilisant les mesures de phase est le calcul preliminaire dit « de 
^initialisation » dans laquelle on resout le probleme des ambiguftes sur le 
nombre de longueurs d'onde. Ce calcul necessite classiquement la 
connaissance prealable d'une position estimee du mobile, qui peut etre 
obtenue en particulier par une methode telle que celle decrite dans les 
35 brevets FR 2 7i5 230 et FR 2 764 708. Cette position estimee est ensuite 



recalee sur la position precise, puis validee au cours de ce calcul 
^initialisation. 

On considere plus particulierement dans la suite I'etape de 
recalage de la position estimee vers une position precise. 
5 La qualite de cette position precise depend notamment de 

I'eloignement du mobile par rapport a la station de reference. 

En effet, les imperfections de la methode differentiate viennent 
d'abord du fait que les signaux radio-satellitaires ne rencontrent pas 
exactement les memes conditions de propagation sur les trajets satellite- 
0 station et satellite-mobile. Les differences des conditions rencontrees, a peu 
pres nulles a proximite immediate de la station, s'accroissent naturellement 
avec I'eloignement. 

Cette difference est principalement due a I'ionosphere qui est 
traversee par les signaux satellite-station et satellite-mobile en des points 
5 differents alors que I'ionosphere n'est pas un milieu homogene. Les mesures 
differentiates basees sur les temps de propagation des signaux satellite- 
station et satellite-mobile sont alors affectees par cette difference. Cette 
difference peut conduire a une erreur de position du mobile par rapport a la 
station de reference pouvant atteindre 1 a plusieurs cm par km 
d'eloignement. Ainsi, pour un eloignement entre la station et le mobile 
superieur a une distance de I'ordre de 10 km, on ne peut garantir la position 
du mobile par rapport a la station de reference avec une precision 
centimetrique. 

Une premiere solution decrite dans le brevet n° 2 764 708 propose 
d'une part de reduire le temps de calcul ^initialisation en particulier le temps 
de calcul d'une position inambigue approchee en utilisant notamment des 
combinaisons lineaires des frequences Remission L1 et L2 des satellites du 
systeme GPS. Elle propose d'autre part de reduire I'erreur ionospherique ; la 
reduction de I'erreur ionospherique s'applique au cours de I'etape de 
recalage. Elle consiste a calculer a partir de la position inambigu# 
approchee, d'une part une position (XL1, YL1, ZL1) pour L1 et d'autre part 
une position (XL2, YL2, ZL2) pour L2, la position precise (X, Y, Z) resultant 
alors de la combinaison lineaire suivante : 
X = (1,65XL1-XL2)/0,65 



Y = (1,65 YL1-YL2)/0,65 
Z=(1,65ZL1-ZL2)/0,65. 

Mais du fait du calcul d'une position sur L1, entache d'une erreur 
ionospherique E et du calcul d'une position sur L2 entache d'une erreur 
5 ionospherique 1,65*E, il n'est pas toujours possible de lever les ambiguites 
lorsque I'erreur ionospherique s'accroit ce qui se produit lorsque la distance 
entre le mobile et la station de reference s'accroit. 

Une autre solution classiquement proposee consiste a mettre en 
place non plus une, mats plusieurs stations de reference constituent ce qui 
10 est communement appele un « Reseau ». Selon cette technique, il est 
possible de connaitre, non seulement les erreurs mesurees en un point 
comme on le fait dans le cadre du « GPS differentiel », mais aussi leur 
gradient devolution dans la zone. On compense done en grande partie I'effet 
des decorellations spatiales des erreurs. Cette solution est efficace, mais elle 
15 est bien sur lourde et couteuse a mettre en csuvre en raison de 
I'infrastructure qu'elle necessite et du cout des communications entre 
stations et mobile. De plus, une telle infrastructure n'existera pas partout. 

II a egalement ete propose une methode basee sur I'exploitation 
du fait que I'erreur ionospherique est fonction de la frequence (en 1/f 2 en 
20 premiere approximation). 

II est alors possible de determiner cette erreur ou de la reduire, 
voire de I'eliminer, en remplacant dans les calculs la frequence f (designee 
par L1 ou L2 dans le cas du systeme GPS) par une combinaison lineaire des 
frequences porteuses des signaux emis par les satellites, e'est-a-dire par une 
25 combinaison lineaire de L1 et L2. 

Le resultat d'une combinaison lineaire de L1, L2 est une nouvelle 
frequence L3 a laquelle correspond une longueur d'onde dite longueur 

d'onde apparente. 

Par exemple, dans le cas du systeme GPS pour lequel L1= 
30 1 ,57542 GHz (ce qui correspond a une longueur d'onde d'environ 19 cm) et 
L2= 1 ,22760 GHz (ce qui correspond a une longueur d'onde d'environ 24 
cm), la combinaison de frequences dite « lono Free », 9L1-7L2 permet 
d'eliminer presque completement I'erreur ionospherique. La longueur d'onde 
apparente correspondante est de 5 cm. 



Cependant cette methods est tres difficile a appliquer du fait de la 
grande difficult a lever les ambiguites sur des longueurs d'onde aussi 
courtes. 
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5 Le but de Invention est done de proposer un precede et un 

dispositif permettant d'obtenir avec une precision centimetrique la position 
d'un mobile par rapport a une station de reference eventuellement eloignee 
de plusieurs dizaines de. kilometres. 

Pour atteindre ce but, I'invention propose un precede de 
10 determination de la position relative d'un mobile par rapport a la position 
connue d'une station de reference, utilisant chacun une antenne de reception 
de signaux radio issus d'une constellation de satellites de positionnement 
emettant sur au moins deux frequences L1 et L2, ce procede comprenant la 
determination periodique pour chacune desdites frequences, d'un jeu de 2p 
pseudo-distances, a savoir p pseudo-distances entre le mobile et les p 
satellites et p pseudo-distances entre la station de reference et les p 
satellites, la fourniture des pseudo-distances a un organe de calcul de 
position, et le calcul par cet organe d'une position relative du mobile par 
rapport a la station de reference a partir d'une part des pseudo-distances et 
d'autre part d'une position estimee Pe du mobile par rapport a la station de 
reference, ce procede etant principalement caracterise en ce que le calcul de 
position relative comprend, pour un jeu donne de 4p pseudo-distances recu 
par i'organe de calcul les etapes suivantes consistant a : 

- a) choisir une combinaison lineaire aL1+bL2 desdites frequences 
25 L1 et L2 parmi une liste predetermine comportant au moins deux 

combinaisons lineaires de frequences, 

- b) calculer les combinaisons lineaires de pseudo-distances 
correspondant a ladite combinaison lineaire et, a partir de ces combinaisons 
lineaires de pseudo-distances et de la position estimee Pe, calculer une 
position relative precise Pp du mobile par rapport a la station de reference, 

- c) choisir dans la liste la combinaison lineaire suivante si elle 
existe, et reiterer dans ce cas I'etape b) en considerant comme position 
estimee, ladite position precise Pp, et en utilisant le meme jeu de 4p pseudo- 
distances, pour obtenir une position relative encore plus precise, 
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- d) reiterer I'etape c) pour toutes les combinaisons lineaires de la 

liste. 

On execute ainsi successivement, a partir du meme jeu de 
mesures de pseudo-distances plusieurs calculs de position du mobile en 
utilisant differentes combinaisons lineaires de frequences, la position estimee 
en entree d'un calcul etant la position calculee a I'etape precedente. 

Dans la technique anterieure un seul calcul etait effectue et la 
combinaison lineaire utilisee dans ce calcul portait sur des positions 
respectivement calculees pour chacune des frequences L1 et L2. 

Une caracteristique importante de I'invention est que les 
combinaisons lineaires de la liste sont determines de telle sorte que les 
longueurs d'onde correspondantes diminuent progressivement et que la 
sensibilite aux erreurs ionospheriques diminue egalement progressivement 
et de facon plus rapide que la longueur d'onde. 

En d'autres mots, la premiere combinaison lineaire de frequences 
est choisie pour que sa longueur d'onde soit grande, afin de faciliter le lever 
d'ambiguite sur le nombre de tours de phase ; mais en contrepartie cette 
premiere combinaison lineaire peut correspondre a une sensibilite importante 
aux erreurs ionospheriques. La position estimee Pe ayant ete amelioree suite 
a ce premier recalage, la deuxieme combinaison lineaire de la liste 
correspond a une longueur d'onde plus faible et a une erreur ionospherique 
encore plus faible. Et ainsi de suite, en utilisant des longueurs d'onde et des 
sensibilites aux erreurs ionospheriques de plus en plus faibles. 

Dans le cadre du systeme GPS, la premiere combinaison peut 
etre la combinaison L1-L2 (a=1, b=-1) ; la derniere peut §tre 9L1-7L2 (a=9, 
b=-7), combinaison connue pour etre pratiquement insensible aux erreurs 
ionospheriques. 

Les combinaisons intermediates sont de preference les suivantes 

(dans I'ordre) : 

3 2L1-L2;3L1-2L2;4L2-3L1. 

Selon un autre aspect caracteristique de I'invention, I'etape b) du 
calcul est effectuee soit en une seule etape utilisant directement les p 
satellites, soit en deux etapes dont la premiere utilise seulement un nombre 
reduit p'' (p'<p) de satellites et la deuxieme utilise les p .satellites. De 

5 deference, c'sst seulement lors de I'utilisation de la premiere combinaison 



lineaire (de longueur d'onde la plus grande) que le calcul est fait en deux 
etapes, les autres combinaisons lineaires etant utilisees selon un calcul en 
une seule etape avec les p satellites. 

Lorsque I'etape b est effectuee en deux etapes, elle comprend 
5 avantageusement les etapes suivantes consistant a : 

b1) calculer une position relative approchee Pa du mobile par 
rapport a la station de reference, a partir de la combinaison lineaire choisie, 
de Pe et dun sous-jeu de 4p' pseudo-distances correspondant a p' satellites 
ou p' est inferieur a p et ou les p' satellites choisis dans la constellation de p 
10 satellites sont ceux qui ont, compte-tenu de la geometrie actuelle de la 
constellation, la plus faible sensibilite a une erreur de la position estimee, 

b2) calculer une position relative precise Pp du mobile par rapport 
a la station de reference a partir de ladite combinaison lineaire, de Pa et du 
jeu complet de 4p pseudo-distances. 
15 Se,on une a^re caracteristique de I'invention, les 2p pseudo- 

distances entre les satellites et la station de reference sont determinees par'i 
la station de reference et envoyees par radio au mobile qui comporte alors. 
des moyens de reception pour recevoir ces pseudo-distances et des 
informations de datation de la mesure de ces pseudo-distances. 
20 Enfin . I'invention a pour objet non seulement le procede de 

determination de position relative d*un mobile par rapport a la position 
connue d'une station de reference dont les grandes lignes viennent d'etre 
decrites, mais aussi un dispositif de determination de la position d'un mobile 
par rapport a une station de reference, capable de mettre en ceuvre ce 
25 procede . Le dispositif selon I'invention comprend au moins, dans le mobile, 
des moyens pour recevoir des signaux de positionnement par satellite, et des 
moyens pour recevoir un jeu de 2p pseudo-distances emises par la station 
de reference et representant les pseudo-distances entre la station de 
reference et p satellites pour au moins deux frequences porteuses differentes 
30 L1 et L2, des moyens de determination periodique d'un jeu de 2p pseudo- 
distances entre le mobile et les p satellites, des moyens de fourniture des 4p 
pseudo-distances a un organe de calcul de position, des moyens de 
stockage d'une liste de combinaisons lineaires des frequences des porteuses 
des signaux de positionnement, des moyens pour effectuer, a partir d'un 
35 meme jeu de 4p pseudo-distances, des calculs successifs de position 



relative du mobile par rapport a la position de la station de reference, chaque 
fois a partir d'une combinaison lineaire differente de frequences choisie dans 
la liste, d'une position estimee Pe et du jeu de 4p pseudo-distances, la 
position estimee lors d'un calcul avec une combinaison lineaire donnee de la 
5 liste etant la position relative calculee a partir de la combinaison lineaire 
precedente de la liste. 

D'autres caracteristiques et avantages de I'invention apparaftront 
a la lecture de la description detaillee qui suit, faite a litre d'exemple non 
1 o iimitatif et en reference aux dessins annexes dans lesquels : 

-la figure 1 represente schematiquement une repartition des 
positions de la station de reference, du mobile et des positions estimee, 
approchee et precise du mobile par rapport a la position de la station de 
reference ; 

15 _|a figure 2 represente schematiquement le dispositif selon 

I'invention ; 

- la figure 3 represente un organigramme des calculs effectues. 

On a represente figure 1, la position R de la station de reference, 
20 eloignee d'une distance D qui peut etre de plusieurs dizaines de kilometres 
d'un mobile dont la position vraie est M. Sont egalement representees sur 
cette figure les positions estimees Pe, approchee Pa et precise Pp du mobile 
auxquelles il sera fait reference plus loin. Le point Pe correspond a la 
position M du mobile estimee a quelques metres pres ( 1 a 2 metres par 

25 exemple), le point Pa a celle du mobile calculee avec une erreur de un a 
quelques decimetres, le point Pp a celle du mobile calculee avec une erreur 
maximum de quelques centimetres. 

La station de reference et le mobile sont respectivement equipes 
d'une anterine de reception 10 et 12 de signaux radio issus des satellites de 

30 positionnement (systeme GPS ou autre) et de moyens de demodulation et de 
traitement des signaux regus. Le mobile calcuie periodiquement des pseudo- 
distances entre sa position et la position des satellites a un instant donne. La 
station calcuie de la meme maniere des pseudo-distances (mesurees au 
meme instant, ou ramenees au meme instant de mesure) entre sa position et 

35 la position des satellites. Comme on I'a deja vu dans le preambule, les 
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pseudo-distances sont classiquement fournies sous forme d'une premiere 
valeur numerique qui definit la position temporelle du code pseudo-aleatoire 
emis par un satellite a un moment donne,. et d'une deuxieme valeur 
numenque qui definit la phase du signal radiofrequence module par ce code 
au meme moment. 

Dans les recepteurs (mobiles ou stations de reference) simples a 
fa.ble precision, un seul signal radiofrequence code est utilise par les circuits 
de reception. Dans les recepteurs plus precis, auxquels la presente invention 
s interesse, au moins deux signaux de frequence porteuse differente L1 et 
L2 en provenance des satellites, sont traites .par le recepteur. Le recepteur 
calcule done periodiquement des pseudo-distances pour chacune des 
frequences L1 et L2. 

II en resulte que pour une mesure de position a un instant donne 
a part.r de p satellites, le calcul de position va faire intervenir un jeu de 4p 
pseudo-distances, done 4p mesures de phase, qui sont respectivement ■ 

aa . ir* • , meSUr6S dS PhaSe 6ntre la Station de reference et p 
satellites, a la frequence L1, 

- p mesures de phase entre le mobile et les p satellites a la 
frequence L1 encore, ' 

- p mesures de phase entre la station de reference et les p 
satellites a la frequence L2, 

- et p mesures de phase entre le mobile et les p satellites a la 

frequence L2. 

La station de reference est equipee de moyens pour emettre vers 
le mobile les valeurs numeriques qu'elle a determines. Le mobile est equipe 
de moyens pour les recevoir (antenne 1 1 sur la figure 2), de sorte que pour 
un calcul de position a un moment donne, le mobile dispose non seulement 
de ses propres mesures de pseudcMiistances mais aussi des mesures de 
pseudo-distances determines par la station de reference. 

Dans le mobile, les signaux provenant de la station de reference 
et les signaux directement recus par le mobile sont utilises par un ensemble 
electronique 14 que comporte le mobile et represents figure 2. On pourrait 
cependant imaginer la solution inverse ou les calculs sont faits par la station 
de reference, le mobile envoyant a celie-ci les pseudo-distances qu'il a 
calculees. Dans ce mode, on connaTt a chaque instant a la station I'endroit 



ou se trouve le mobile, ce qui peut etre interessant dans certaines 
applications. 

Cet ensemble electronique 14 comprend tout d'abord un circuit de 
reception de positionnement par satellite 16 capable de determiner des 
5 mesures de pseudo-distances entre son antenne et plusieurs satellites. 

Le circuit de reception 16 foumit periodiquement (par exemple 
toutes les 100 millisecondes) un jeu de 2p pseudo-distances lorsque p 
satellites sont en vue directe avec Tantenne 12. II les transmet a un organe 
de calcul de position relative, 18, qui recoil par ailleurs le jeu de 2p pseudo- 
10 distances correspondant aux 2p pseudo-distances recues des meme p 
satellites par Tantenne 10 de la station de reference et transmises au mobile 

parl'antenne 11. 

L'organe de calcul de position relative 18 a pour fonction de 
determiner la position precise Pp de Tantenne 12 du mobile par rapport a 
1 5 Tantenne 1 0 de la station de reference. 

L'organe de calcul 18 est programme pour effectuer les calculs 
desires et est relie aux organes peripheriques necessaires selon les 
applications : afficheur 20, clavier 22, iecteur/enregistreur de donnees ou de 
programmes 24, moyen de transmission filalre ou radio vers un utilisateur, ou 
simple interface de sortie vers une ligne de transmission du resultat des 
calculs. 

Les moyens de calcul prevus dans l'organe de calcul sont aptes a 
effectuer les operations qu'on va detainer ci-apres, a partir d'un jeu de 4p 
pseudo-distances recues a un instant donne par le mobile. 

Le calcul de position est differentiel, c'est-a-dire qu'on determine 
la position du mobile par rapport a la station de reference (que le calcul soit 
effectue par le mobile ou par la station). On peut done considerer dans la 
suite que le calcul de position consiste a calculer la position du mobile a 
partir de mesures differentielles de pseudo-distances, en faisant des 
30 differences entre pseudo-distances mesurees a la station de reference et au 
mobile. Le calcul fait intervenir une notion de doubles differences de pseudo- 
distances (differences entre couples de satellites) que Ton precisera plus 
loin. 

Globalement, en utilisant la notion classique de double difference, 
35 le principe de calcul differentiel est le suivant : 
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- on calcule les positions des satellites a ('instant de mesure t 
grace aux ephemerides des satellites ; 

- on determine, a partir des mesures differentielles de pseudo- 
distances, les distances differentielles entre la station de reference et Ie 

5 mobile selon les axes de visee des satellites. On obtient des distances qu'on 
peut considerer globalement comme les projections, Ie long de ces axes de 
la distance D entre Ie mobile et la station de reference ; ce sont des 
distances mesurees ; 

- on calcule pour Ie meme instant de mesure les distances Ie long 
io des memes axes entre la station de reference et une position estimee Pe du 

mobile ; ce sont des distances estimees ; 

- on determine selon chaque axe la difference entre la distance 
mesuree et la distance estimee, qu'on peut appeler grandeur d'ecart ou 
« innovation » selon cet axe ; 

15 - on calcule a partir de ces grandeurs d'ecart, a I'aide de la matrice 

des cosmus directeurs representant les directions des axes de visee des 

estTmee 8 ' ** ^ * P ° Slti ° n m ° bi ' e 6t ,a position 

- on rajoute a la position estimee initiate les ecarts calcules et on 
obtient une position calculee du mobile qui est soit une position definitive Pp 
soit une nouvelle position estimee Pe en vue d'une etape ulterieure de calcui 
comme on Ie verra plus loin. 

Comme on I'expliquera plus loin, les doubles differences 
considerees dans la presente invention sont elaborees, non pas a partir des 
frequences directement mesurees, mais a partir de. combinaisons lineaires 
des pseudo-distances mesurees a la frequence L1 et des pseudo-distances 
mesurees a la frequence L2. .. 

Un element essentiel de la presente invention reside dans Ie fait 
qu on execute successivement, a partir du meme jeu de 4p mesures de 
pseudo-distances plusieurs calculs de position du mobile en utilisant 
differentes combinaisons lineaires des phases correspondent aux 
frequences L1 et L2. Cela revient a calculer des phases fictives d'une 
frequence porteuse qui serait fictivement la combinaison lineaire des 
frequences L1 et L2. Plusieurs combinaisons lineaires sont successivement 
utilisees pour Ie meme jeu de 4p pseudc-distances ; Ie calcul utilise les 
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pseudo-distances et une position estimee et il aboutit a une position calculee 
pour une combinaison lineaire donnee ; la position calculee pour cette 
combinaison Hneaire sert de position estimee pour un calcul suivant utilisant 
une autre combinaison lineaire. La succession de combinaisons lineaires, 

5 correspondant a differentes longueurs d'onde apparentes, est telle que la 
precision du calcul augmente progressivement. Les combinaisons lineaires 
de frequences sont choisies d'un calcul a Tautre notamment de maniere a ce 
que ia demi-longueur d'onde apparente correspondante soit superieure ou 
egale a I'erreur de position resultant du calcul precedent, tout en diminuant 

i o I'erreur ionospherique. 

Cornme on le verra plus loin, on utilise dans les calculs de 
position, des doubles differences etablies pour des couples de satellites 
determines de fagon a augmenter encore la longueur d'onde apparente en 
prenant en compte la geometrie de la constellation des satellites (la longueur 

1 5 d'onde apparente d'un couple de satellites est d'autant plus grande qu'il est 
vu sous un angle etroit). 

En combinant judicieusement le choix des combinaisons lineaires 
de frequences et des couples de satellites, on peut ainsi passer d'une 
position estimee de precision decimetrique, voire metrique, a la position de 

20 precision centimetrique recherchee. Les combinaisons lineaires sont choisies 
dans une liste dont un exemple est donne plus loin. Cette liste est stockee 
par exemple dans I'organe 24 pour etre utilisee par Torgane de calcul 18. 

Dans le cas du GPS t les combinaisons lineaires s'appliquent aux 
frequences L1 et L2. Le resultat d'une combinaison lineaire de L1, L2 est 

25 une nouvelle frequence L3 a laquelle correspond une longueur d'onde dite 
longueur d'onde apparente. Le circuit de reception 16 fournit a I'organe de 
calcul 18, la phase cpu du signal emis par le satellite pour L1 et la phase cpi_2 
du signal emis par le satellite pour L2, ces deux phases etant representatives 
de la distance entre le satellite et le mobile ; I'organe de calcul 18 pourra 

30 alors utiliser non pas la phase cpu ou la phase cpL2 mais une phase <pi_3 
correspondant a la combinaison lineaire suivante 

L3 = a.L1 + b.L2 selon la formule cpu = a. <pu + b. <pL2 
De meme, la station de reference envoie au mobile des phases 
representatives, pour chaque satellite, de la distance entre le satellite et la 



station de reference et la meme combinaison lineaire peut etre utilisee pour 
determiner une phase apparente a la longueur d'onde L3. 

Les nombres a et b sont les coefficients de la combinaison lineaire 

choisie. 

5 . Pour chaque combinaison lineaire, on sait calculer non seulement 

la longueur d'onde apparente (done la distance au-dela de laquelle la mesure 
de phase devient ambigue) mais aussi la sensibilite aux erreurs de 
propagation ionospherique. Si on prend comme valeur d'erreur 
ionospherique de reference I'erreur existant a la frequence L1, on peut 

10 calculer la valeur correspondante de I'erreur pour chaque combinaison 
lineaire. Si par exemple I'erreur sur L1 est de 1 cm/km d'eloignement, pour un 
eloignement de 10 km a plusieurs dizaines de kilometres, I'erreur sur L1 
atteindrait 10 cm a plusieurs dizaines de cm. 

Pour corriger cette erreur, on pourrait imaginer d'utiliser 

15 directement une combinaison lineaire de frequences donnant une tres faible 
erreur de propagation ionospherique ; la combinaison dite « tono-Free » 9L1- 
7L2 a une tres faible erreur ; mais elle correspond a une longueur d'onde 
apparente tres courte (5cm) de sorte qu'il est impossible de lever son 
ambiguTte de phase a partir d'une position estimee eloignee. 

20 Pour le GPS et a titre d'exemple prefere, on a etabli la liste 

ordonnee suivante de combinaisons lineaires, avec la longueur d'onde 
apparente correspondante et le coefficient d'erreur ionospherique 
(relativement a la valeur de reference unitaire pour la frequence L1). Les 
longueurs d'onde apparentes sont decroissantes, et les coefficients d'erreur 

25 ionospherique decroissent encore plus rapidement 

L1-L2 connue de I'homme du metier sous le nom de « Wide-Lane », 
qui correspond a une longueur d'onde apparente d'environ 86 cm avec un 
ratio de I'erreur ionospherique par rapport a la frequence L1 de base de 1,3 ; 
2L1-L2, qui correspond a une longueur d'onde apparente d'environ 16 

30 cm avec un ratio de I'erreur ionospherique par rapport a la frequence L1 de 
base de 0,56 ; 

3L1-2L2, qui correspond a une longueur d'onde apparente d'environ 
13 cm avec un ratio de I'erreur ionospherique par rapport a la frequence L1 
de base de 0,3 ; 
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4L2-3L1, qui correspond a une longueur d'onde apparente d'environ 
11 cm avec un ratio de l'erreur ionospherique par rapport a la frequence L1 
de base de 0,1 ; 

9L1-7L2, denommee comme on I'a deja vu combinaison « lono-Free » 
qui correspond a une longueur d'onde apparente d'environ 5 cm avec un 
ratio de l'erreur ionospherique par rapport a la frequence L1 de base proche 
dezero. 

Bien sOr cette sequence de combinaisons lineaires nest quun 
exemple adapte au GPS dans sa situation actuelle : elle differera notamment 
dans le cas de frequences L1 et L2 differentes de celles citees (futures 
evolutions du GPS, autres systemes radio-satellitaires tels que Galileo) ; elle 
peut aussi ne contenir que certaines de ces combinaisons lineaires et/ou en 

incorporer d'autres. 

Le principe est qu'avec un jeu de 4p pseudo-distances combmees 
lineairement, on peut effectuer un calcul en utilisant une position estimee et 
la premiere combinaison lineaire de la liste, avec un faible risque d'amb.gu.te 
de phase du fait de la longueur d'onde importante de cette combinaison ; on 
aboutit a une premiere position calculee plus precise que la position est.mee. 
Cette position calculee sert de position estimee pour un autre calcul fart avec 
les memes 4p pseudo-distances mais combinees selon la combinaison 
suivante de la liste, correspondent a une longueur d'onde apparente plus 
faible mais encore suffisamment grande pour ne pas introduce d'ambigu.te 
de phase compte-tenu de l'erreur ionospherique subsistant suite au prem.er 
calcul On aboutit a une nouvelle position calculee toujours plus precise. 
Cette position plus precise sert de position estimee pour un troisieme calcul 
avec le meme jeu de pseudo-distances combinees selon la tro.s.eme 
combinaison lineaire de la liste, et ainsi de suite. A chaque etape la position 
est amelioree du fait de la reduction de l'erreur ionospherique, et 
Amelioration de cette position permet de choisir, dans I'etape suivante et 
sans risque d'erreur d'ambiguite, une longueur d'onde plus courte qui avec 
une reduction toujours plus importante de l'erreur ionospherique permet de 
rapprocher encore la position calculee de la position vraie. 

La figure 3 represente schematiquement les grandes etapes des 

calculs effectues. 
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Un jeu de 4p pseudo-distances est fourni a I'organe de calcul 
(etape 1) qui a determine par ailleurs une position estimee Pe et qui choisit 
la premiere combinaison lineaire de frequences de ia liste predeterminee 
(§tape2). 

5 Facultativement, le calcul d'une position a partir d'une position 

estimee et d'une combinaison lineaire donnee peut se faire en deux grandes 
etapes: 

-une position approchee Pa est calculee a partir de la position 
estimee Pe, de la combinaison lineaire de frequences choisie et d'un jeu de 
10 4p' pseudo-distances (etape 3a). 

- une position precise Pp est ensuite calculee a partir du point 
approche Pa, de la meme combinaison lineaire de frequences choisie et du 
jeu complet de 4p pseudo-distances (etape 3b). 

Dans une variante de I'invention, on peut calculer directement en 
une etape la position precise Pp a partir de la position estimee Pe, de la 
combinaison lineaire de frequences choisie et d'un jeu complet de 4p 
pseudo-distances ; on peut par exemple retenir le calcul en deux etapes pour 
la premiere combinaison lineaire et en une etape pour les combinaisons 
lineaires suivantes. 

Si toutes les combinaisons lineaires de la liste n'ont pas encore 
ete utilisees (etape 4), on choisit la combinaison lineaire suivante de la liste 
et on reitere les etapes 3a et 3b en considerant comme position estimee la 
position Pp calculee resultant de I'etape 3b (etape 5). 

Si toutes les combinaisons lineaires de la liste ont ete utilisees on 
25 valide (etape 6) la position Pp resultant du dernier calcul en determinant 
classiquement la valeur de coherence de la position calculee, pour rejeter les 
solutions qui ne respecteraient pas des criteres de coherence minimaux. A 
Tissue de cette etape de validation, la position Pp est la position du mobile 
etablie avec une precision centimetrique. 

Dans le cas ou la position Pp n'est pas validee, on execute les 
etapes decrites a partir d'un autre jeu de 4p pseudo-distances. 

On va donner maintenant, a titre d'exemple, le detail du calcul des 
etapes 3a et 3b qi/on peut faire pour I'obtention d'une position precise Pp a 
partir d'une position estimee Pe. 
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La phase d'obtention d'une position approchee Pa (etape 3a) 
utilise une position estimee initiate Pe ; et elle utilise avantageusement un 
sous-jeu de 4p' pseudo-distances choisies dans le jeu de 4p pseudo- 
distances. 

5 On effectue le calcul des positions par un traitement en double 

difference, sur les 4p' pseudo-distances. 

Le traitement dit « en double difference » consiste a travailler non 
pas directement a partir des differences pour deux satellites entre pseudo- 
distances mais a partir de differences entre la station de reference et le 
1 o mobile, de la difference pour deux satellites entre pseudo-distances. 

Ces couples de satellites sont choisis en fonction de leur 
sensibilite aux erreurs de position, cette sensibilite dependant de la 
geometrie entre les satellites et le point de mesure. La geometrie de visibility 
des p satellites a I'instant de la mesure etant connue grace aux ephemerides 
15 et a la position estimee, on sait classer les couples de satellites dans I'ordre 
de leur sensibilite croissante aux erreurs de position. On ne prend que les p' 
satellites correspondant aux couples les moins sensibles, ce qui permet 
d'augmenter la longueur d'onde apparente. Plus precisement, on calcule les 
doubles differences du type 
20 DDij = (Dim-Djm)-(Dir-Djr), dans lesquelles : 

Dim est la pseudo-distance du mobile au satellite de rang i 
Djm est la pseudo-distance du mobile au satellite de rang j 
Dir est la pseudo-distance de la station de reference au satellite 

de rang i 

25 Djr est la pseudo-distance de la station de reference au satellite 

de rang j 

Les differences du type Dim-Djm ou Dir-Djr permettent d'eliminer 
les erreurs communes aux satellites ( differences d'horloge entre les 
satellites et le recepteur). 
30 Les differences entre ces differences, ou doubles differences DDij, 

permettent d'eliminer les erreurs dues a la propagation atmospherique ou 
ionospherique. 

Les doubles differences sont calculees d'une part pour la 
frequence L1 (difference DDijl), d'autre part pour la frequence L2 (difference 
35 DDij2). On calcule siors, pour une premiere combinaison lineaire de 



frequences, L3 = aL1+bl_2, une combinaison lineaire de doubles differences 
CLij - aDDijl +bDDij2. Ce sont ces combinaisons lineaires qui vont §tre 
utilisees, et non les classiques doubles differences, pour le calcul de 
position ; elles represented une phase a la longueur d'onde apparente de la 
5 frequence Active L3 et sont exprimees ici en distances. 

Les combinaisons de doubles differences CLij sont comparees a 
des combinaisons similaires, calculees et non mesurees, a partir de la 
position estimee initiate Pe. La difference qui resulte de cette comparaison 
est appelee INNOVij, representant I'ecart entre I'estimation et la mesure. 
10 Ces ecarts sont relies aux ecarts de longitude, latitude, et altitude 

DL, DG et DA entre la position estimee (ici la position initiate Pe) et la 
position calculee (ici la position calculee approchee Pa) par des equations 
du type 

INNOVij = DL[cos(Evi)cos(Azi)-cos(Evj)cos(Azj)] 
1 5 + DG[cos(Evi)sin(Azi)-cos(Evj)sin(Azj)J 

+DA[sin(Evi)-sin(Evj)] 

ou Evi, Evj sont les elevations des satellites i et j, et Azi, Azj leurs - 

azimuts. 

Un calcul simple, ou un calcul matriciel avec minimisation des 
20 erreurs par la technique des moindres carres si on a plus de 4 satellites, < 
permet de determiner DL, DG, DA qui represented des ecarts entre position 
mesuree et position estimee. Ces ecarts sont ajoutes a la longitude, la 
latitude, et I'altitude de la position estimee Pe pour obtenir une position 
approchee Pa. 

25 A P artir de cette position approchee Pa on effectue une deuxieme 

etape de calcul (etape 3b). La deuxieme etape est tres semblable a la 
premiere, mais 

- elle utilise toutes les 4p pseudo-distances, c'est-a-dire celles 
correspondant a tous les couples de satellites, 

30 - ©'I© utilise comme position estimee non pas la position 

estimee initiate Pe mais la position approchee Pa 

- elle utilise pour le calcul d'une position precise Pp un calcul 
matriciel avec un nombre d'equations en general superieur au nombre 
d'inconnues (le nombre p de satellites etant suppose superieur a 4) ; la 

35 determination des ecarts DL, DG, DA entre position estimee et position 



calculee peut se fairs alors classiquement par une methode des moindres 
carres (la position determinee par calcul est celle qui minimise la valeur 
quadratique moyenne des residus). 

La position Pp est ensuite utilisee dans un nouveau calcul comme 
position estimee, a la place de la position Pe, et en faisant intervenir la 
combinaison lineaire suivante de la liste. 

Progressivement, par iteration jusqu'au bout de la liste, on arrive a 
une valeur de plus en plus precise et non ambigue de la position du mobile. 
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REVINDICATIONS 

1 . Procede de determination de la position relative d'un mobile par 
rapport a la position connue d'une station de reference, utilisant chacun une 
antenne de reception de signaux radio issus d'une constellation de satellites 
de positionnement emettant sur au moins deux frequences L1 et L2, ce 

5 procede comprenant la determination periodique pour chacune desdites 
frequences, d'un jeu de 2p pseudo-distances, a savoir p pseudo-distances 
entre le mobile et les p satellites et p pseudo-distances entre la station de 
reference et les p satellites, la fourniture des pseudo-distances a un organe 
de calcul de position, et le calcul par cet organe d'une position relative du 

10 mobile par rapport a la station de reference a partir d'une part des pseudo- 
distances et d'autre part d'une position estimee Pe du mobile par rapport a la 
station de reference, ce procede etant principalement caracterise en ce que 
le calcul de position relative comprend, pour un jeu donne de 4p pseudo- 
distances regu par I'organe de calcul, les etapes suivantes consistant a : 

15 - a) choisir une combinaison lineaire aL1+bL2 desdites frequences 

L1 et L2 parmi une liste predeterminee comportant au moins deux 
combinaisons lineaires de frequences, 

- b) calculer les combinaisons lineaires de pseudo-distances 
correspondant a ladite combinaison lineaire et, a partir de ces combinaisons 

20 lineaires de pseudo-distances et de la position estimee Pe, calculer une 
position relative precise Pp du mobile par rapport a la station de reference, 

- c) choisir dans la liste la combinaison lineaire suivante si elle 
existe, et reiterer dans ce cas i'etape b) en considerant comme position 
estimee, ladite position precise Pp, et en utilisant le meme jeu de 4p pseudo- 

25 distances, pour obtenir une position relative encore plus precise, 

- d) reiterer I'etape c) pour toutes les combinaisons lineaires de la 

liste. 

2. Procede selon la revendication precedente, caracterise en ce 
30 que les combinaisons lineaires de la liste sont determinees de telle sorte 

que, d'un calcul a I'autre, les longueurs d'onde correspondantes diminuent 
progressivement et que la sensibilite aux erreurs ionospheriques diminue 
egalement progressivement. 



ler depot 




3. Procede selon Tune quelconque des revendications 
precedentes, caracterise en ce que la premiere combinaison lineaire de la 
liste est la combinaison L1-L2 (a=l, b=-1)et/ou la derniere combinaison 
lineaire de la liste est la combinaison 9L1-7L2 (a=9, b=-7), L1 et L2 etant les 

5 frequences d'emission des satellites du systeme GPS. 

4. Procede selon la revendication precedente, caracterise en ce 
que les combinaisons lineaires intermediates sont de preference les 
suivantes (dans I'ordre) : 

10 2L1-L2 (a=2, b=-1); 3L1-2L2 (a=3, b=-2) ; 4L2-3L1 (a=4, b=-3). 

5. Procede selon Tune quelconque des revendications 
precedentes, caracterise en ce que I'etape b) comprend les deux etapes 
suivantes consistant a : 

15 b1) calculer une position relative approchee Pa du mobile par 

rapport a la station de reference, £ partir de la combinaison lineaire choisie, 
de Pe et d'un sous-jeu de 4p' pseudo-distances correspondent a p 1 satellites 
ou p' est inferieur a p et ou les p 1 satellites choisis dans la constellation de p 
satellites sont ceux qui ont, compte-tenu de la geometrie actuelle de la 

20 constellation, la plus faible sensibilite & une erreur de position estimee, 

b2) calculer une position relative precise Pp du mobile par rapport 
a la station de reference a partir de ladite combinaison lineaire, de Pa et du 
jeu complet de 4p pseudo-distances. 

25 6. Proced§ selon la revendication precedente, caracterise en ce 

que les etapes b1) et b2) ne sont realisees que pour la premiere 
combinaison lineaire de la liste, une seule etape faisant intervenir les 4p 
pseudo-distances etant realisee pour les autres combinaisons lineaires de la 
liste. 

30 

7. Procede selon Tune quelconque des revendications 
precedentes, caracterise en ce que les 2p pseudo-distances entre les 
satellites et la station de reference sont determinees par la station de 
reference et envoyees par radio au mobile qui comporte alors des moyens 



de reception pour recevoir ces pseudo-distances et des informations de 
datation de la mesure de ces pseudo-distances. 

8. Dispositif de determination de la position d'un mobile par 

5 rapport a une station de reference, comprenant au moins, dans le mobile, 
des moyens pour recevoir des signaux de positionnement par satellite et des 
moyens pour recevoir un jeu de 2p pseudo-distances emises par la station 
de reference et representant les pseudo-distances entre la station de 
reference et p satellites pour au moins deux frequences porteuses differentes 

10 L1 et L2, des moyens de determination periodique d'un jeu de 2p pseudo- 
distances entre le mobile et les p satellites, des moyens de fourniture des 4p 
pseudo-distances a un organe de calcul (18) de position, des moyens de 
stockage d'une liste de combinaisons lineaires des frequences des porteuses 
des signaux de positionnement, des moyens pour effectuer, a partir d'un 

15 meme jeu de 4p pseudo-distances, des calculs successifs de position 
relative du mobile par rapport a la position de la station de reference, chaqije 
fois a partir d'une combinaison lineaire differente de frequences choisie dans 
la liste, d'une position estimee Pe et du jeu de 4p pseudo-distances, (a 
position estimee lors d'un calcul avec une combinaison lineaire donnee de la 

20 liste etant la position relative calculee a partir de la combinaison lineaire 
precedente de la liste. 
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Rua 


THALES Intellectual Property J 
13, Avenue du President Salvador Allende | 




Code postal etville 


94117 lARCUEIL 


Society d*appartenance (facultatif) 
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Prenoms 




Adresse 


Rue 






Coda postal etville 
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Nom 





Prenoms 



Coda postal et ville 



Societe d'appartenance (facultatif) 
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